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Abstract 

We study the evolution of fronts in the Klein-Gordon equation when the nonlin- 
ear term is non-homogeneous. Extending previous works on homogeneous nonlinear 
terms, we describe the derivation of an equation governing the front motion, which is 
strongly nonlinear, and, for the two-dimensional case, generalizes the damped Born- 
Infeld equation. We study the motion of one- and two-dimensional fronts, finding that 
the dynamics is richer than in the homogeneous reaction term case. 



1 Introduction 

In the last few years, partial differential equations with discrete nonlinearities have been 
used to model phenomena in fields ranging from physics to biology, including the study of 
*E-mail: horacio@cs.brandeis.edu 
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pinning of dislocation motions in crystals, breathers in nonlinear crystal lattices, Josephson 
junction arrays and the biophysical description of calcium release waves [?, ?, ?, ?, ?, ?, ?, 
?, ?, ?, ?, ?, ?, ?, ?, ?, ?, ?]. Recently, the discrete one- dimensional stationary version of 
the Klein-Gordon equation 

(ptt + 7 <Pt = D (j) xx + a Y,6(x- x k ) [f{4>) + h], (1) 

k 

with 7 = has been analyzed by Flach and Kladko [?] (see also references therein). In (|l|) 
<f> is an order parameter, the non-negative constant 7 is the dissipation coefficient and the 
positive constants D and a are the diffusion coefficient and the amplitude of the discrete 
nonlinearity, respectively. The function / is a bistable function (the derivative of a double 
well potential having the two equal minima); i.e., a real odd function with positive maximum 
equal to <f>*, negative minimum equal to —(f)* and precisely three zeros in the closed interval 
[a_,a+] located at a_, ao and a+. For simplicity and without lost of generality we will 
consider in our analysis a_ = — 1, a = and a + = 1. The prototype example is f(<f>) = 
(<f) — <f) 3 )/2. The constant h, assumed to be small in absolute value, specifies the difference of 
the potential minima of the system; i.e., /(</>) + h is the derivative of a double well potential 
with one local minimum and one global minimum. Note that ([I]) reduces to the Klein-Gordon 
equation when J2k $( x ~ x k) is replaced by a constant with appropriate rescaling. In [?] a 
first order perturbation calculation for the heteroclinic orbits of the corresponding stationary 
kink solution for ([[]) was presented. Kink solutions, connecting the two local minima of the 
double well potential, were also obtained for the sine-Gordon case, f((f>) = — sin(</>), and the 
Klein-Gordon case, f(<f>) = (0 — <f> 3 )/2. Both are particular cases of the function f((f>) as 
defined above. Note that the sine- Gordon case is equivalent to the derivative of a double 
well potential in a restricted domain of definition. 

In this manuscript, we study the dynamics of kinks for a quasi-discrete version of the Klein- 
Gordon equation 
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4> tt + 1 4>t = D A<p + aP(x,y) [f{4>) + h], (2) 

in a bounded region Q C R n , n = 1, 2 with smooth boundary <9fi for Neumann boundary 
conditions on dQ. Equation (^) reduces to ([I]) when (3 is one-dimensional and j3(x) = 
J2k$i x ~ x k)- Altough the analysis presented below will be valid for a general class of 
positive different iable functions (3, we have in mind some particular cases which approximate 
a distribution of discrete nonlinearities for large r], a positive constant defined below. 

Case 1) There is a sequence of points on the real line, Xk, k = 1, . . . , N, with iV finite or 
infinite, where the function (3 reaches a maximum, 

N 

p{x) = e-^ x - Xk)2 . 

k=l 

Case 2) There is a sequence of lines in the plane, yk, k — 1, . 
where the function j3, independent of x, reaches a maximum, 

N 

P(x,y) = Yl e-^-y^ 2 . (4) 

k=l 

Case 3) There is a sequence of points in the plane, (x^, yj), k — 1, . . . , N, j — 1, . . . , M with 
N and M finite or infinite, where the function (3 reaches a maximum, 

N M 

P(x,y) = EE e(x-Xk, V - Vj] v)> where a(x, y; rj) = e~ v{x +y K (5) 

k=ij=i 

Case 4) There is a sequence of circles in the plane, p — p^, k — 1, . . . , N, with iV finite or 
infinite, and where p represents the radial polar coordinate, where the function (3 reaches a 
maximum, 

N 

M = E e-^) 2 . (6) 
k=i 



(3) 

N, with iV finite or infinite, 
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We refer to the points Xk and (xk, Vj),k = 1, . . . , N, j = 1, . . . , M as quasi-discrete (QD) sites 
and to the stripes y = Vk and circles p = pk k = 1, . . . , N as quasi-semi-discrete (QS) sites. 
We define d to be the minimum distance between two adjacent QD or QS sites. Note that 
the function (3 can be chosen to depend not only on the spatial variable but also on t. The 
specific form of /3(x, y, t) will depend on the particular model. One might, for example, have 
the product of a spatially dependent function (3(x,y) with a probabilistic time-dependent 
function. 

For (0) we define the following dimensionless variables and parameters 



x= d' y =d> t = — (7) 



< D 1 7 d , 2 c- h 

7 = — ■= rj = rj d , h 



and 



a d ^JD ' ' e 

Substituting (0) and @ into (0) and dropping the " from the variables and parameters 

e 2 <f> tt + e 2 7 t = e 2 A0 + y) [/(0) + e (9) 

We will consider the case 0<e<l; i.e., when diffusion is slow, d is large or a is large, and 
there is a small dissipation. 

The homogeneous version of @, 

e 2 0« + e 2 7 0t = e 2 A0 + /(0) + h, (10) 

possesses a travelling kink solution. The point on the line (for n = 1) or the set of points 
in the plane (for n — 2) for which the order parameter <p vanishes are called the interface 
or front of the system. For (TO) the front moves according to an extended version of the 
Born-Infeld equation [?, ?] 



(1 - s 2 ) s xx + 2s x s t s xt - (1 + si) s tt - 7 s t (1 + s 2 x - s 2 ) - h (1 + s\ - s t 2 )§ = 0, (11) 

where h, proportional to h, will be denned later. Planar fronts moving according to ( |TTD 
with 7 = h — (no dissipation and both phases with equal potential) move with a constant 
velocity equal to their initial velocity. For other values of 7 or h, fronts move with a velocity 
that asymptotically approaches —h/{^ 2 + h 2 )^ as long as the initial velocity is bounded from 
above by 1 in absolute value. Linear perturbations to these planar fronts decay, either in a 
monotonic or an oscillatory way, to zero as t — > 00 [?] . Circular interfaces moving according 
to ( PH) with h > shrink to a point in finite time [?, ?]. If h < 0, then circles shrink to 
points for some initial conditions and for others they grow unboundedly. Neu [?] showed 
that for 7 = h = 0, closed kinks can be stabilized against collapse by the appearance of short 
wavelength, small amplitude waves. For the more general case, two situations are possible. 
Either linear perturbations to a circle decay and curves shrink to a point in finite time or 
they are still present at the shrinkage point of the circle. Note that Equation ([□]) expressed 
in terms of its kinematic and geometric properties reads [?] 

< ^ + ^v(l-v 2 )-K(l-v 2 ) + h(l- v 2 )^ = 0, (12) 

where k is the curvature of the front and dv/dt is the "Lagrangian" time derivative of v which 
is calculated along the trajectory fo the interfacial point moving with the normal velocity v 

[?]■ 

One of the goals of this papers is to determine whether the dynamic behavior of kinks in (§) 
differs from its homogeneous (discrete) nonlinearity counterpart (|i~0|). For the overdamped 
version of equation (Pp, which is a parabolic bistable equation, it has been shown that there 
are essential differences between the homogeneous and non-homogeneous (discrete) cases, in 
that latter exhibits propagation failure [?, ?, ?, ?, ?]. 
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In Section 2 we present an equation of motion for the front in equation ([5]), and we describe 
briefly the method by which it was derived. This equation generalizes equation (|TTD with 
the strong nonlinearity accounting for the influence of the function (3 on the front motion. 
In Section 3 we study the evolution of one-dimensional fronts. We show that for h = the 
function (3 acts as a "potential function" for the motion of the front; i.e., a front initially 
placed between two maxima of (3 asymptotically approaches the intervening minimum. When 
h 7^ 0, fronts that start between two maxima of j3 asymptotically approach an equilibrium 
point determined by h and /3, producing a kink propagation failure. In Section 4 we study 
the evolution of two-dimensional fronts with radial symmetry. We find that when there is 
no dissipation circles can shrink to a point in finite time, grow unboundedly or their radius 
can oscillate, depending on the initial conditions. When dissipation effects are present, 
the oscillations decay spirally or not depending of the value of 7. The final result is the 
stabilization of the circular domain of one phase inside the other phase. Our conclusions 
appear in Section 5. 

2 Front Dynamics: The Equation of Motion 

For (§) and the law of motion of the interface in two dimensions is given by 
(1 - s 2 ) s xx + 2 s x s t s xt - (1 + s 2 x ) s tt - 7 s t (1 + s 2 x - s 2 )- 

_ p y (x, s) - & (a s) g + (3 t (x, s)s t (1 + ^ _ ^ _ h g) (1 + s 2_ s 2 ) | = 0; (13) 

where y = s(x,t) is the Cartesian description of the interface and h is proportional to h 
as will be explained later. Equation flT3] ) was obtained by carrying out a non-rigorous but 
self-consistent singular perturbation analysis for e < 1, treating the interface as a moving 
internal layer of width 0(e). We focused on the dynamics of the fully developed layer, and 
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not on the process by which it was generated. The method that we applied is similar to that 
used in [?] and [?] for the study of the evolution of kinks in both the nonlinear wave equation 
([10|) and the Allen-Cahn equation with quasi-discrete sources of reaction (the overdamped 
version of (§)). The basic assumptions made were: 

- For small e > and all t G [0, T], the domain Q can be divided into two open regions £l+(t; e) 
and Q-(t,e) by a curve T(t;e), which does not intersect dQ. This interface, defined by 
r(t; e) := {x E Q : <j>(x,t; e) = 0}, is assumed to be smooth, which implies that its curvature 
and its velocity are bounded independently of e. 

- There exists a solution <ft(x, £; e) of @, defined for small e, for all x G Q and for all t G [0, T] 
with an internal layer. As e — > this solution is assumed to vary continuously through the 
interface, taking the value 1 when x G Q + (t;e), —1 when x G fi_(£, e), and varying rapidly 
but smoothly through the interface. 

- The curvature of of the front is small compared to its width. 

As a first stage in the derivation of equation (|13|) we define near the interface a new variable 
z = (y — s)/e which is 0(1) as e — > oo and then express equation (|9]) in terms of this new 
variable. Next we expand <p and (3 asymptotically in a power series in e and substitute these 
expansions into the differential equation. After equating the coefficients of corresponding 
powers of e, we obtain two equations. The first can be reduced to an equation of the 
type + /($°) = which has to satisfy $°(0) = and $°(±1) = ±1, giving a kink 
solution. Here $° represents the leading order term of the order parameter in terms 
of z. The second problem is a linear non-homogeneous second order ODE. Equation ( |13|) 
is obtain by applying the solvability condition (Fredholm alternative) after defining h : = 
h [$°(+oo) - $°(-oo)] / JZzi^zfdz. Note that for /(0) = (0 - 3 )/2 (Ginsburg-Landau 
theory), <3>°(z) = tanh | and h = 3 h whereas for /(</>) = sin0 (sine- Gordon), $>°(z) = 
4 tan~ 1 e 2: — n and h = | h. 
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3 Front Motion in ID 



For a one-dimensional system, equation fllTf ) reads 



s tt + 7 s t (1 - # + ^r(l -sl) + h j3l(a) (1 - s 2 )§ = 0. (14) 

We concentrate on functions (3 of the form (|3p, altough the same analysis can be done for a 
general differentiable function. We define u = s and v = s t obtaining 

v* = v 1 

(15) 

v t = - 7 v (1 - t; 2 ) - ^gy(l -v 2 )-h fi{u) (1 - t; 2 )!. 

The fixed points of (p~5|) are (tt , 0), where wo satisfies #(w) = /?'(«) + 2 h fti (u) = 0. 
The trace, r, and determinant, A, of matrix of the linearized system are r = —7 and 
A = [(3"(u ) (3(u )-(3' 2 (u )} I 2 (3 2 (u ) + h f3'(u )/[2 P*(uq)], respectively. If h = then the 
fixed points are the maxima (unstable) and minima (stable) of (3{u). Thus, a front initially 
between two maxima of (3 will move and asymptotically approach the intervening minimum. 
When there is no dissipation, this behavior is in contrast with the homogeneous case ( O ) 
where, as was pointed out in the introduction, fronts move with a constant velocity equal 
to their initial velocity. In the non-homogeneous case, we can predict the final position of 
the front from the structure of (3. In order to understand the behavior of g(u) as h increases 
above zero we consider a function (3{u) with a single peak at 0; i.e., (3{u) = e~ vu2 . This 
function will approximate the more general (§) if rj ^> 1, so that the influence of peaks on 
one another is very small. In this case g{u) = —2e~ r]u2 [ri u — h e~ n ~?~\. For h — 0, g(u) 
vanishes at u = and it is positive for u > and negative for u < 0. As A moves away from 
zero, u, the root of g(u), will be given by the solution of 77 u — h e~ R 2~ = 0, an equation that 
always has a solution. If h > 0, then x > 0, and g(u) is positive for x > x and negative for 
x < x. If h < 0, then x < 0. As an illustration, We can see the shape of g(u) as h increases 
in Figure |l|. In summary, as h increases or decreases the behavior of the front is similar to 
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the case h = 0, in contrast to the classical homogeneous case ( |TT| ) where, as noted in Section 
1, fronts with an initial velocity whose absolute value is bounded from above by 1, move 
with a velocity that asymptotically approaches —h/{^ 2 + h 2 )^. 



4 Front Motion in 2 D 

The analysis of front motion in two dimensions governed by (|T3|) with a function (3 of type 
(d) reduces to the analysis of one- dimensional front motion, and we shall not consider this 
case further. For radially symmetric functions, (3 = /3(p), and radially symmetric fronts, 
equation (|l3|) for the radial coordinate p of the front reads 

Pu + (7 Pt + ~ p ) (1 - P\) + (1 - P 2 t ) + PHp) h (1 - pDl (16) 

We define u = p and v = p t obtaining 



(17) 



u t = v, 

vt = -b v + I + j^jgj + h p\ (u) (1 - v 2 )^(l - v 2 ). 
The lines v = ±1 are trajectories of (|H]) in the corresponding phase plane. They define 
a region D with the property that every curve starting in this region remains inside it for 
all future time. We confine our analysis to u > 0. We analyze here the case h = 0. The 
fixed points of ([[7]) are (w , 0), where u are solutions (if they exist) of 2 (3{u) + u (3'{u) = 0. 
The trace, r, and the Determinant, A, of the matrix of the linearized system are given by 
T = - 7 and A = -l/u 2 a + [f3"(u ) (3(u ) - (3' 2 (u )} / 2 f3 2 (u ) respectively. The simplest 
case is (3{p) = e~ v( - p ~ Pl ^ 2 for a given p\ > 0. For this case uq = [r] p\ + \j ffp\ + 4 r\ j /2 77 
and (wo,0) is a saddle point. For (|J) with 77 = 50 and iV = 2, pi = 0.5 and P2 = 1.5, 
we calculated the fixed points of (|TTD using the Newton-Raphson method with a tolerance 
of 0.0001. They are z x = 0.537228, z 2 = 0.999165 and z 3 = 1.513217. The corresponding 
values of A are A(^) = -53.464832, A(z 2 ) = 1196.824463 and A(z 3 ) = -50.436714. Then 
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Z\ and z% are saddle points, and z 2 is stable. Since the discriminant, A = r 2 — 4 A of z 2 , is 
—4787.297852 < 0, z 2 is a neutrally stable center for 7 = 0, a stable spiral for < 7 < 70 
and a stable node for 7 > 70, where 70 is the value of 7 for which A = 0. 

For the case 7 = 0, dividing the second equation in (|T7|) by the first and solving one obtains 

c 2 u 2 (3(u) + v 2 = 1, (18) 

where c 2 = (1 — v 2 )/{u 2 f3{ui)) and (ui, vi) are the initial conditions. In Figure we present a 
graph of flT8| ) for (3 given by (f|) with N = 2, p 1 = 0.5, p 2 = 1.5 77 = 50 (a) and 77 = 10 (b). We 
observe that there are shrinking trajectories, oscillatory trajectories and growing trajectories 
in contrast with the homogeneous case where all trajectories are shrinking trajectories given 
by c 2 u 2 + v 2 = 1 with c 2 = (1 — v 2 )/u 2 [?]. As 77 decreases, the oscillatory trajectories 
dissapear, leaving growing trajectories, which will ultimately vanish as rj — > 0. 

In order to study more generally the behavior of the system away from the fixed points we can 
see in Figure |3] the phase plane for 7 = (a) and 7 = 1 (b) . The dashed lines are the nullclines 
of the system and the "o" are its fixed points. The trajectories were calculated solving 
(|T7|) using a Runge-Kutta method of order four. We observe there that there are different 
situations according to the initial conditions. In Figure Q-a (no dissipation), trajectories A, 
B and C correspond to circles that shrink to a point in finite time. If their initial velocity 
is positive, then their radius grows initially to a value bounded by z\ before shrinkage takes 
place. Trajectories D and J also correspond to a circles that finally shrink to a point in finite 
time. In the case of D, although the initial conditions are close to those of trajectory C, 
the dynamics is very different. In addition to shrinkage, trajectories can display unbounded 
growth, represented by trajectory G, and periodic behaviour, represented by trajectories E 
and F. Trajectories H and I also correspond to circles growing unboundedly, but if the 
initial velocity is negative they shrink to a valued bounded from below by z 3 and then they 
start growing. In Figure |]-b (7 = 1) we see that trajectories A, B and C correspond to circles 
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that shrink to points in finite time after growing to a radius bounded by z\. Trajectory D 
also shrinks to a point in finite time, but it grows initially to a radius bounded from above 
by Z3 and from below by zi- As we pointed out before, as a consequence of dissipation 
(7 7^ 0) Z2 is a stable spiral. We see that trajectory E spirals into z 2 , and there are no 
longer periodic trajectories. For N > 2 we expect the phase plane analysis to be similar to 
that presented here. In contrast with the homogeneous nonlinear wave equation, where any 
circular front shrinks to a point in finite time, the nonhomogeneous version (|T3|) presents a 
very rich dynamics with periodic motion and stabilization of circular domains of one phase 
inside the other. 

In the absense of dissipation there are two "forces" responsible for the motion of the front: the 
curvature of the circular front, 1/p, and the "potential function" (3. For initial conditions 
near enough to the minimum of (3 the two "forces" balance and oscillations are possible. 
When dissipation is present that "balance" is lost, and the oscillations decay. 

5 Conclusions 

In this manuscript we have presented equation (|T3|) as governing the evolution of a fully 
developed front in a nonhomogeneous version of the nonlinear wave equation, when e <C 
1. This equation generalizes the damped version of the Born-Infeld equation ([□]) to include 
the effects of stronger nonlinearities and accounts for the influence of the nonhomogeneous 
nonlinear term on the motion of the front. The motion of interfaces according to ( pT3[ ) 
is qualitatively different from that of the homogeneous counterpart given by ([11]). This 
difference arises primarily from the fact that the function (3 acts as a "potential function" 
for the motion of the front. For the one dimensional case, an initial front initially placed 
between two maxima of (3 (which for a homogeneous nonlinear term will move with a velocity 
that asymptotically approaches —hji^j 1 + h 2 )^ as long as the initial velocity bounded from 
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above by 1 in absolute value asymptotically approaches a point depending on h and on the 
structure of (3. For the radially symmetric two-dimensional case, the dynamics is richer than 
in the homogeneous counterpart, where for h — circles shrink to point in finite time. In 
the absence of dissipation, circles can shrink to a point in finite time, grow unboundedly 
or their radius oscillates, depending on the initial conditions. When dissipation effects are 
present, the oscillations decay, spirally or not, depending on the value of 7. The final result 
is the stabilization of a circular domain of one phase inside the other phase. 

The evolution of circular interfaces in more complicated arrangement of QD sites and the 
evolution of more complicated fronts like convex closed curves calls for further research. We 
hope to address this questions in a forthcoming paper. 

Acknowledgement: We thank the Chemistry Division of the National Science Foundation 
for support of this work. 
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Captions 
Figure [I]: 

a) Graph of (3{u) for 77 = 1000, X\ = 2, x 2 = 1, x 3 = and x 4 = —1, x 5 = —2. 

b) Graph of g(u) for 77 = 1000, h — 0, x\ — 2, X2 — 1, #3 = and £4 = — 1, £5 = —2. 

c) Graph of g(u) for 77 = 1000, /i = 10, xi = 2, X2 = 1, x% — and £4 = —1, X5 = —2. 

d) Graph of g(u) for 77 = 1000, h = 20, 37 = 2, x 2 = 1, x 3 = and x 4 = —1, x 5 = —2. 
The points x^ are the maxima of f3(u). 

Figure ||: 

Graph of ((TJ) with /3 given by (g) with jV = 2, pi = 0.5, p 2 = 1.5 and 

a) 77 = 50. 

b) 7/ = 10. 

Figure |3|: 

Phase plane for fllTD with /3 given by @ with N = 2, pi = 0.5, p 2 = 1.5, 77 = 50 and h = 0. 
The dashed lines are the nullclines of the system and the "o" are its fixed points. 

a) 7 = 0. 

b) 7 = l. 
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